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A detection or nondetection of primordial non-Gaussianity by using the cosmic microwave back- 
ground Radiation (CMB) data is crucial not only to discriminate inflationary models but also to 
test alternative scenarios. Non-Gaussianity offers, therefore, a powerful probe of the physics of the 
primordial universe. The extraction of primordial non-Gaussianity is a difficult enterprise since 
several effects of non-primordial nature can produce non-Gaussianity. Given the far-reaching con- 
sequences of such a non-Gaussianity for our understanding of the physics of the early universe, it is 
important to employ a range of different statistical tools to quantify and/or constrain its amount 
in order to have information that may be helpful for identifying its causes. Moreover, different 
indicators can in principle provide information about distinct forms of non-Gaussianity that can 
be present in CMB data. Most of the Gaussianity analyses of CMB data have been performed by 
using part-sky frequency, where the masks are used to deal with the galactic diffuse foreground 
emission. However, full-sky map seems to be potentially more appropriate to test for Gaussianity 
of the CMB data. On the other hand, masks can induce bias in some non-Gaussianity analyses. 
Here we use two recent large-angle non-Gaussianity indicators, based on skewness and kurtosis of 
large-angle patches of CMB maps, to examine the question of non-Gaussianity in the available full- 
sky five-year and seven-year Wilkinson Microwave Anisotropy Probe (WMAP) maps. We show that 
these full-sky foreground-reduced maps present a significant deviation from Gaussianity of different 
levels, which vary with the foreground-reducing procedures. We also make a Gaussianity analysis 
of the foreground-reduced five-year and seven-year WMAP maps with a KQ75 mask, and compare 
with the similar analysis performed with the corresponding full-sky foreground-reduced maps. This 
comparison shows a significant reduction in the levels of non-Gaussianity when the mask is em- 
ployed, which provides indications on the suitability of the foreground-reduced maps as Gaussian 
reconstructions of the full-sky CMB. 

PACS numbers: 98.80.Es, 98.70.Vc, 98.80.-k 



I. INTRODUCTION 

A key prediction of a number of simple single-field 
slow-roll inflationary models is that they cannot gener- 
ate detectable non-Gaussianity of the cosmic microwave 
background (CMB) temperature fluctuations within the 
level of accuracy of the Wilkinson Microwave Anisotropy 
Probe (WMAP) Q. There are, however, several infla- 
tionary models that can generate non-Gaussianity at a 
level detectable by the WMAP. These non-Gaussian sce- 
narios comprise models based upon a wide range of mech- 
anisms, including special features of the inflation poten- 
tial and violation of one of the following four conditions: 
single fleld, slow roll, canonical kinetic energy, and initial 
Bunch-Davies vacuum state. Thus, although convincing 
detection of a fairly large primordial non-Gaussianity in 
the CMB data would not rule out all inflationary mod- 
els, it would exclude the entire class of stationary mod- 
els that satisfy simultaneously these four conditions (see, 
e.g., Refs. 0t3)- Moreover, a null detection of deviation 
from Gaussianity would rule out alternative models of the 
early universe (see, for example, Refs. Q). Thus, a de- 
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tection or nondetection of primordial non-Gaussianity in 
the CMB data is crucial not only to discriminate (or even 
exclude classes of) inflationary models but also to test al- 
ternative scenarios, offering therefore a window into the 
physics of the primordial universe. 

However, there are various non-primordial effects that 
can also produce non-Gaussianity such as, e.g., unsub- 
tracted foreground contamination, unconsidered point 
sources emission and systematic errors d-Q. Thus, the 
extraction of a possible primordial non-Gaussianity is 
not a simple endeavor. In view of this, a great deal of 
effort has recently gone into verifying the existence of 
non-Gaussianity by employing several statistical estima- 
tors (for related articles see, e.g., Refs. 0). Differ- 
ent indicators can in principle provide information about 
multiple forms of non-Gaussianity that may be present 
in WMAP data. It is therefore important to test CMB 
data for deviations from Gaussianity by using a range 
of different statistical tools to quantify or constrain the 
amount of any non-Gaussian signals in the data, and ex- 
tract information on their possible origins. 

A number of recent analyses of CMB data performed 
with different statistical tools have provided indications 
of either consistency or deviation from Gaussianity in the 
CMB temperature fluctuations (see, e.g., Ref. @). In a 
recent paper jTl| we proposed two new large-angle non- 
Gaussianity indicators, based on skewness and kurtosis 
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of large-angle patches of CMB maps, which provide mea- 
sures of the departure from Gaussianity on large angular 
scales. We used these indicators to search for the large- 
angle deviation from Gaussianity in the three and five- 
year single frequency maps with a KQ75 mask, and found 
that while the deviation for the Q, V, and W masked 
maps are within the 95% expected values of Monte-Carlo 
(MC) statistically Gaussian CMB maps, there is a strong 
indication of deviation from Gaussianity (;» 95% off the 
MC) in the K and Ka masked maps. 

Most of the Gaussianity analyses with WMAP data 
have been carried out by using CMB temperature fiuc- 
tuation maps (raw and clean) in the frequency bands Q, 
V and W or some combination of these maps. In these 
analyses, in order to deal with the diffuse galactic fore- 
ground emission, masks such as, for example, KQ75 and 
KpO have been used. 

However, sky cuts themselves can potentially induce 
bias in Gaussianity analyses, and on the other hand full- 
sky maps seem more appropriate to test for Gaussianity 
in the CMB data. Thus, a pertinent question that arises 
is how the analysis of Gaussianity made in Ref. llj is 
modified if whole-sky foreground-reduced CMB maps are 
used. Our primary objective in this paper is to address 
this question by extending the analysis of Ref. [ll| in 



three different ways. First, we use the same statistical in- 
dicators to carry out a new analysis of Gaussianity of the 
available full-sky foregro und-reduced five-year and seven- 
year CMB maps [12|-[l5|. Second, since in these maps the 
foreground is reduced through different procedures each 
of the resulting maps should be tested for Gaussianity. 
Thus, we make a quantitative analysis of the effects of 
distinct cleaning processes in the deviation from Gaus- 
sianity, quantifying the level of non-Gaussianity for each 
foreground reduction method. Third, we study quan- 
titatively the consequences for the Gaussianity analysis 
of masking the foreground-reduced maps with the KQ75 
mask. An interesting outcome is that this mask lowers 
significantly the level of deviation from Gaussianity even 
in the foreground-reduced maps, rendering therefore in- 
formation about the suitability of the foreground-reduced 
maps as Gaussian reconstructions of the full-sky CMB. 



II. NON-GAUSSIANITY INDICATORS 

The chief idea behind our construction of the non- 
Gaussianity indicators is that a simple way of accessing 
the deviation from Gaussianity distribution of the CMB 
temperature fluctuations is by calculating the skewness 
S — ^:i/a^, and the kurtosis K = ^i^/a^ — 3 from the 
fluctuations data, where /is and /i4 are the third and 
fourth central moments of the distribution, and a is its 
variance. Clearly calculating S and K from the whole 
sky temperature fluctuations data would simply yield two 
dimensionless numbers, which are rough measures of de- 
viation from Gaussianity of the temperature fluctuation 
distribution. 



However, one can go further and obtain a great number 
of values associated to directional information of devia- 
tion from Gaussianity if instead one takes a discrete set 
of points {j = 1, . . . , A'c} homogeneously distributed on 
the celestial sphere S*^ as the center of spherical caps of a 
given aperture 7 and calculate Sj and Kj from the CMB 
temperature fluctuations of each spherical cap. The val- 
ues Sj and Kj can then be taken as measures of the 
non-Gaussianity in the direction {9j^ (pj) of the center of 
the spherical cap j . Such calculations for the individual 
caps thus provide quantitative information (2Nc values) 
about possible violation of Gaussianity in the CMB data. 

This procedure is a constructive way of deflning 
two discrete functions S and K (deflned on S*^) from 
the temperature fluctuations data, and can be formal- 
ized through the following steps (for more details, see 
Ref. [ni): 

i. Take a discrete set of points {j = 1, . . . , iVd homo- 
geneously distributed on the CMB celestial sphere 
as the centers of spherical caps of a given aper- 
ture 7; 

ii. Calculate for each spherical cap j the skewness (Sj) 
and kurtosis (Kj) given, respectively, by 



(1) 



and 
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where A^p is the number of pixels in the j cap, 
Ti is the temperature at the i**^ pixel, Tj is the 
CMB mean temperature in the j^^ cap, and a is 
the standard deviation. Clearly, the values Sj and 
Kj obtained in this way for each cap can be viewed 
as a measure of non-Gaussianity in the direction of 
the center of the cap {6j ,(j>j); 

iii. Patching together the Sj and Kj values for each 
spherical cap, one obtains our indicators, i.e., dis- 
crete functions 5* = S{9, (j)) and K — K{9, </>) de- 
fined over the celestial sphere, which can be used to 
measure the deviation from Gaussianity as a func- 
tion of the angular coordinates {9,<j)). The Moll- 
weid projection of skewness and kurtosis functions 
S = S{9, (f>) and K = K{9, (/>) are nothing but skew- 
ness and kurtosis maps, hereafter we shall refer to 
them as S*— map and A"— map, respectively. 

Now, since S = S{9, (p) and K = K{9, (j)) are functions 
defined on 5^ they can be expanded into their spherical 
harmonics in order to have their power spectra Si and Ki^. 
Thus, for example, for the skewness indicator S — S{9, cj)) 
one has 



(3) 



e=0 ra=- 



3 



and can calculate the corresponding angular power spec- 
trum 
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which can be used to quantify the angular scale of the 
deviation from Gaussianity, and also to calculate the sta- 
tistical significance of such deviation. Obviously, similar 
expressions hold for the kurtosis K — K{6, (j)). 

In the next section we shall use the statistical indica- 
tors S = S{9, (j)) and K = K{9, (j)) to test for Gaussianity 
the available foreground-reduced maps obtained from the 
five-year WMAP data. 



III. NON-GAUSSIANITY 

A. Foregound-reduced maps 

The WMAP team has released high angular resolution 
five-year maps of the CMB temperature fluctuations in 
the five frequency bands K (22.8 GHz), Ka (33.0 GHz), 
Q (40.7 GHz), V (60.8 GHz), and W (93.5 GHz). They 
have also produced a full-sky foreground- reduced Internal 
Linear Combination (ILC) map which is formed from a 
weighted linear combination of these five frequency band 
maps in which the weights are chosen in order to mini- 
mize the galactic foreground contribution. 

It is well known that the first-year ILC map is inap- 
propriate for CMB scientific studies However, in the 
five-year (also in the three-year and seven-year) version 
of this map a bias correction has been implemented as 
part of the foreground cleaning process, and the WMAP 
team suggested that this map is suitable for use in large 
angular scales (low t) analyses although they admittedly 
have not performed non-Gaussian tests on this version 
of the ILC map [H, 113 • Notwithstanding the many 
merits of the five-year ILC procedure, some cleaning fea- 
tures of this ILC approach have been considered, and 
two variants have been proposed recently. In the first 
approach the frequency dependent weights were deter- 
mined in harmonic space [13], while in the second the 
foreground is reduced by using needlets as the basis of the 
cleaning process [l3|. Thus, two new full-sky foreground- 
cleaned maps have been produced with the WMAP five- 
year data, namely the harmonic ILC (HILC) |13| an d the 
needlet ILC (NILC) (for more details see Refs. [il,[il). 

In the next section, we use the full-sky foreground- 
reduced ILC, HILC and NILC maps with the same 
smoothed 1° resolution (which is the resolution of the 
ILC map) as the input maps from which we calculate the 
S — S{d, (f)) and K ~ K{9, 0) maps, and then we com- 
pute the associated power spectra in order to carry out 
a statistical analysis to quantify the levels of deviation 



from Gaussianity. ""^ 



B. Analysis and results 

In order to minimize the statistical noise, in the cal- 
culations of skewness and kurtosis maps (S*— map and 
ii'— map) from the foreground-reduced maps, we have 
scanned the celestial sphere with spherical caps of aper- 
ture 7 = 90°, centered at 12 288 points homogeneously 
generated on the two-sphere by using the HEALPix 
code [ll]. In other words, the point-centers of the 
spherical caps are the center of the pixels of a homo- 
geneous pixelization of the S"^ generated by HEALPix 
with A^sidc = 32. We emphasize, however, that this pix- 
elization is only a practical way of choosing the centers 
of the caps homogeneously distributed on S^. It is not 
related to the pixelization of the above-mentioned ILC, 
HILC and NILC input maps that we have utilized to cal- 
culate both the S and K maps from which we compute 
the associated power spectra. 

Figures [T] and [D show examples of S and K maps ob- 
tained from the foreground-reduced NILC full-sky and 
KQ75 maps. The panels of these figures clearly show re- 
gions with higher and lower values ('hot' and 'cold' spots) 
of S{9,(j)) and K{9,(l)), which suggest large-angle multi- 
pole components of non-Gaussianity. We have also cal- 
culated similar maps (with and without the KQ75 mask) 
from the ILC and HILC maps. However, since these maps 
provide only qualitative information, to avoid repetition 
we only depict the maps of Figs. [T] and [5] merely for illus- 
trative purpose. 

In order to obtain quantitative information about the 
large angular scale (low i) distributions for the non- 
Gaussianity S and K maps obtained from the available 
full-sky foreground-reduced five-year maps, we have cal- 
culated the (low £) power spectra Si and Kg for these 
maps. The statistical significance of these power spec- 
tra is estimated by comparing with the corresponding 
multipole values of the averaged power spectra Si and 
Kg calculated from maps obtained by averaging over 
1 000 Monte-Carlo-generated statistically Gaussian CMB 
maps.^ Throughout the paper the mean quantities are 
denoted by overline. 

Before proceeding to a statistical analysis, let us de- 
scribe with some detail our calculations. For the sake 
of brevity, we focus on the skewness indicator S, but 
a completely similar procedure was used for the kur- 



The ILC, HILC and NILC maps are available for download from: 
http: / /lambda. gsfc.nasa.gov/product / map/ dr3 / ilc_map_get.cfm] 
http:/ /www.nbi.dk/f^jkim/hilc/ and 



http: / / www.apc.univ-paris7.fr / APC_CS/Recherche/ Adamis /cmb_wmap-en.p] 
Each Montc-Carlo scrambled map is a stochastic realization of 
the WMAP best-fitting angular power spectrum of the ACDM 
model, obtained by randomizing the temperature components 
a»™ within the cosmic variance limits. 
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FIG. 1: Skewness indicator maps calculated from the five-year foreground-reduced NILC full-sky (left panel) and KQ75 masked 
(right panel) maps. 
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FIG. 2: Kurtosis indicator maps calculated from the five-year foreground-reduced NILC full-sky (left panel) and KQ75 masked 
(right panel) maps. 



tosis indicator K. We generated 1 000 MC Gaussian 
{scrambled) CMB maps, which are then used to generate 
1 000 skevifness S*— maps, from which we calculate 1 000 
power spectra: {S^} (i — 1, ••• ,1000 is an enumera- 
tion index, and = 1, • • • , 10 ). In this way, for each 
fixed multipole component Sg^f^^^^ we have 1 000 mul- 
tipole values from which we calculate the mean value 
Se = (1/1000) SI . From this MC process we have 

at the end ten mean multipole values Si, each of which 
are then used for a comparison with the corresponding 
multipole values Se (obtained from the input map) in 
order to evaluate the statistical significance of the mul- 
tipole components Se- To make this comparison easier, 
instead of using the angular power spectra Si and Ki 
themselves, we employed the differential power spectra 
\Se — Se\ and \Ke — Ke\, which measure the deviation 
of the skewness and kurtosis multipole values (calculated 
from the foreground-reduced maps) from the mean mul- 
tipoles Si and Kg (calculated from the Gaussian maps). 
Thus, for example, to study the statistical significance of 
the quadrupole component of the skewness from HILC 
map S^^^'^ (say) we calculate the deviation \S^^^'^ -^2], 
where the mean quadrupole value S2 is calculated from 



the i = 1, • • • , 1 000 quadrupole values of the MC Gaus- 
sian maps. 

Figure [3] shows the differential power spectra calcu- 
lated from full-sky five-year foreground-reduced maps, 
i.e., it displays the absolute value of the deviations from 
the mean angular power spectrum of the skewness Si 
(left panel) and kurtosis Kg (right panel) indicators 
for £ — 1, ••• ,10, which is a range of multipole val- 
ues needed to investigate the large-scale angular char- 
acteristics of the S and K maps. This figure shows a 
first indication of deviation from Gaussianity in five-year 
foreground-reduced ILC, HILC and NILC maps in that 
the deviations \Si — Stl and \Ki — Ki\ for these maps are 
not within 95% of the mean MC value. 

To obtain additional quantitative information regard- 
ing the deviation from Gaussianity, we can also calcu- 
late the percentage of the deviations \S\ — Si\ calculated 
from 1 000 MC Gaussian maps, which are smaller than 
\Si — Si\ obtained from each foreground-reduced map. 
This calculations are made in detail in the Appendix VK\ 
Thus, for example, we have for the full-sky NILC, HILC 
and ILC maps, respectively, that ~ 99.999%, - 99.999%, 
and 99.900% of the multipole values S'g obtained from the 
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FIG. 3: Differential power spectrum of skewness \Se — Si\ (left) and kurtosis \Kt — Kt\ (right) indicators calculated from 
the full-sky foreground-reduced ILC, BlILC, and NILC maps obtained from the WMAP five-year data. The 68% and 95% 
confidence levels are indicated, respectively, by the dashed and dash-dotted lines. 





for 5, 


X" for Ke 


HILC 


4 625 


301 665 


ILC 


35.7 


2 368 


NILC 


7.1 


160.3 



TABLE I: Results of the test to determine the goodness of fit for 5*^ and Ki multipole values calculated from the full-sky 
foreground-reduced HILC, ILC, and NILC maps as compared to the expected multipoles values from the Caussian MC maps. 



MC maps are closer to the mean £"5 than the value 6*5 cal- 
culated from the data, i.e. from each of the foreground- 
reduced maps. This indicates how unlikely are the oc- 
currences of the values obtained from these foreground- 
reduced maps for the multipole in the set of values 
of 5*5 from MC simulated maps. In other words, the 
probability of occurrence of the S5 values (in the set of 
MC values) for the NILC, HILC and ILC maps is only 
O(10-3)%, 0(10-3)% and 0(10-^)%, respectively. Sim- 
ilarly, the probability of occurrence of K2, for example, is 
0(10-'^)% for all these foreground-reduced maps, while 
for K5 are respectively O(10-i)% (NILC), 0(10-^)% 
(HILC) and O{10-^)% (ILC). In Tables |IV] and IVJ of the 
Appendix 1X1 we collect together the probability of occur- 
rence of each of the values Si and {£ = 1, ••• ,10) 
calculated from S and K maps obtained from the full- 
sky NILC, HILC and ILC maps. In Tables El and IVHl 
we present these probabilities calculated from the same 
input maps but now with KQ75 mask."^ The compar- 



^ We emphasize that, throughout this paper, in the implementa- 
tion of the mask we do not take T = for the temperature 
fluctuation of the pixels inside the masked region. This would 
clearly induce non-Gaussian contribution. In our scan of the 
CMB sky when the spherical cap move into the masked area the 
pixels of the cap inside the masked do not contribute to the val- 
ues of the indicators in the center of the cap. In these cases, the 
values Sj and Kj for a j cap are calculated with small number 



ison of Table |TVl with Table El and of Table EI with 
Table En makes apparent the role of the KQ75 mask in 
reducing the level of deviation from Gaussianity (see the 
Appendix |X1 for more details). 

Although the set of 'local' (fixed £) estimates collected 
together in the tables of Appendix |^ gives an indication 
of deviation from Gaussianity as measured by each mul- 
tipole component to have an overall assessment of low 
£ power spectra Sg and calculated from each CMB 
foreground-reduced map, we have performed a test to 
find out the goodness of fit for Si and Kg multipole val- 
ues as compared to the expected multipole values from 
the MC Gaussian maps. In this way, we can obtain 
one number for each foreground-reduced map that collec- 
tively ('globally') quantifies the deviation from Gaussian- 
ity. For the power spectra Se and Kg we found that the 
values given in Table U for the ratio x^/dof (dof stands 
for degrees of freedom) for the power spectra calculated 
from HILC, ILC and NILC full-sky input maps. Clearly 
a good fit occurs when x^/dof ^ 1. Moreover, greater 
are the x^/dof values, the smaller the x^ probabilities, 
that is the probability that the multipole values Se and 
Kg and the expected MC multipole values agree. Thus, 
regarding the skewness indicator Table U shows that the 
HILC presents the greatest level of deviation from Gaus- 
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FIG. 4: DifTerential power spectrum of skewness \Si — Se\ (left) and kurtosis \Kt — Ke\ (right) indicators calculated from the 
five-year foreground-reduced KQ75 masked ILC, IflLC and NILC maps. The 68% and 95% confidence levels are indicated, 
respectively, by the dashed and dash-dotted lines. 





for St [KQ75] 


for Ki [KQ7.5] 


HILC 


4.7 


4.2 


ILC 


1.2 


0.4 


NILC 


1.4 


1.1 



TABLE II: Results of the x^ test to determine the goodness of fit for S£ and Ki multipole values calculated from the foreground- 
reduced HILC, ILC, and NILC maps with a KQ75 mask as compared to the expected multipoles values from the Gaussian MC 
masked maps. 





for Si [KpO] 


X^ for Ke [KpO] 


HILC 


58.7 


101.9 


ILC 


1.9 


6.5 


NILC 


4.5 


17.9 



TABLE III: Results of the test to determine the goodness of fit for Se and Ki multipole values calculated from the S and 
K maps obtained from the foreground-reduced HILC, ILC, and NILC maps with a KpO mask as compared to the expected 
multipoles values from the Gaussian MC masked maps. 



sianity {x^/doi ^ 1), as captured by the indicator S, 
while the NILC map has the lowest level. 

Regarding the deviation from Gaussianity as detected 
by the kurtosis indicator K, Table U shows again that the 
HILC presents the largest deviation followed by the ILC 
and NILC. To the extent that x^/dof is considerably 
greater than one, all these full-sky foreground-reduced 
maps also present a significant deviation from Gaussian- 
ity as captured here by the kurtosis indicator. 

The above results of our statistical analysis given in 
Fig. [3] and gathered together in Table U (and also sup- 
ported by Tables IIVI and |V] of Appendix |X| show a sig- 
nificant deviation from Gaussianity in five-year full-sky 
foreground-reduced (ILC, NILC and HILC) maps as de- 
tected by both the skewness and the kurtosis indicators 
S and K. A pertinent question that arises here is how 
this analysis of Gaussianity for the full-sky foreground- 
reduced maps is modified if one uses the KQ75 mask. 



which was recommended by the WMAP team for tests 
of Gaussianity of the five-year band maps. Furthermore, 
the combination of the full-sky and mask analyses should 
provide information on the reliability of the foreground- 
reduced maps as appropriate reconstructions of the full- 
sky CMB. 



Figure m shows the power spectra \S£ — Se\ (left) and 
\Ki — Ki\ (right) calculated from five-year foreground- 
reduced KQ75 masked maps. This figure along with 
Fig. [3] show a significant reduction in the level of de- 
viation from Gaussianity when the foreground-reduced 
ILC, HILC, and NILC maps are masked. To quantify 
this reduction we have recalculated x^/dof for these in- 
put maps with the KQ75 mask, and have collected the 
results in Table [Hi The comparison of Table U and Ta- 
ble [H] shows quantitatively the reduction of the level of 
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Gaussianity for the case of CMB masked maps.^ 

In the above analyses we have followed the five-year 
WMAP recommendation for tests of Gaussianity and 
thus used the mask the KQ75, which is slightly more 
conservative than the KpO {theKQ75 sky cut is 28.4% 
while the KpO cut is 24.5%). A pertinent question at 
this point is how the above results are modified if the 
less conservative KpO mask is used. We have examined 
this issue by calculating the power spectra Si and Kg and 
the x^/dof from S and K maps obtained from the ILC, 
NILC, and HILC input maps with the KpO mask. The 
result of this analysis is given in Table IIIIl 

A comparison between Tables [TT] and IIIIl shows that in 
general the value of x^/dof increases for both indicators 
when the less conservative mask KpO is used. We note 
that the changes in x^/dof values are greater for the 
HILC, though. 

The comparison between Fig.[3]and Fig.|4j and TablesH] 
and Table |lT] along with the tables of the Appendix |X] 
clearly provides quantitative information on the suitabil- 
ity of the foreground-reduced maps as Gaussian recon- 
structions of the full-sky CMB, and makes apparent the 
relevant role of the mask KQ75 in reducing significantly 
the level of non- Gaussianity in these foreground-reduced 
maps. 

The calculations of our non- Gaussianity indicators re- 
quire the specification of some quantities whose choice 
could in principle affect the outcome of our calculations. 
To test the robustness of our scheme, hence of our results, 
we studied the effects of changing in the parameters em- 
ployed in the calculation of our indicators. We found 
that the S and K angular power spectra do not change 
appreciably as we change the resolution of CMB temper- 
ature maps used and the number of point-centers of the 
caps with values 768, 3 072 and 12 288 (see Ref. ^ for 
more details on the robustness of this method) . 

Concerning the robustness of the above analyses with 
the KQ75 mask some additional words of clarification 
are in order here. First, we note that the calculations of 
the S*— maps and ii'— maps by scanning the CMB masked 
maps sometimes include caps whose center is within or 
close to the KQ75 masked region. In these cases, the 
calculations of the S and K indicators are made with a 
smaller number of pixels, which clearly introduce addi- 
tional statistical noise as compared to the full-sky map 
cases. In order to minimize this effect we have scanned 
the CMB masked sky with spherical caps of aperture 
7 = 90°, and for the sake of uniformity we have used 
caps with the same aperture for the full-sky maps. We 
note, however, that full-sky foreground-reduced analysis 
does not change significantly if one uses smaller apertures 
as, for example, 7 ~ 60°. 



* Incidentally, this reduction is also revealed through (and agrees 
with) the c ompa rison of Table |IV] with Table jV] and of Table Ivll 
with Table rVlTl 



IV. CONCLUDING REMARKS 

The detection or nondetection of primordial non- 
Gaussianity in the CMB data is essential to discriminate 
or even exclude classes of infiationary models. It can also 
be used to test alternative scenarios of the primordial 
universe. There are, however, several non-primordial ef- 
fects that can also produce non-Gaussianity. This makes 
the extraction of a possible primordial non-Gaussianity 
a rather difficult endeavor. Since different indicators can 
in principle provide information about distinct forms of 
non-Gaussianity, it is important to test CMB data for 
non-Gaussianity by using different estimators to quantify 
and/or constrain its amount in order to extract informa- 
tion about their possible sources. 

Most of the Gaussianity analyses of CMB data have 
been performed with frequency band maps. In these 
studies, to deal with the galactic diffuse foreground emis- 
sion, masks have been employed. However, a full-sky 
foreground-reduced map seems to be potentially more 
appropriate to test for Gaussianity the CMB data.^ The 
five-year version of the ILC map has been suggested as 
a full-sky map suitable for large angular scales analy- 
ses , even though the WMAP team has not performed 
a battery of non-Gaussianity tests on this map (I^ . 

In this paper we have performed an analysis of Gaus- 
sianity of the available five-year full-sky foreground- 
reduced maps. To this end, we have used two new non- 
Gaussianity indicators based on skewness and kurtosis 
of large-angle patches of CMB maps, which provide a 
measure of departure from Gaussianity on large angu- 
lar scales [ll|. We have shown that the full-sky five- 
year foreground-reduced maps (ILC, HILC and NILC) 
present a significant deviation from Gaussianity, which 
varies with the foreground-reducing procedures. We have 
established which of these full-sky foreground-reduced 
maps exhibit the highest and the lowest level of non- 
Gaussianity. 

We have also masked the foreground-reduced maps 
with KQ75 and KpO masks and performed a quantita- 
tive analysis of deviation from Gaussianity of these maps. 
The comparison of the full-sky and masked analyses (see 
Fig. H and Fig. H and Tables U III and |lll| shows a sig- 
nificant reduction in the levels of non-Gaussianity when 
the masks are employed, which in turn provides indica- 
tions on the suitability of the foreground-reduced maps 
as Gaussian reconstructions of the full-sky CMB. 

Finally, when we were in the process of rewriting a re- 
vised version of this paper, by taking into account the 
referee's recommendations, the seven-year WMAP CMB 



^ In reality, the full-sky map seems to be the most suitable for a 
number of other issues, including the test of statistical isotropy, 
the search for evidence of a North-South asymmetry in CMB 
data, and signatures of a possible nontrivial cosmic topology, for 
example. 
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I NILC [full-sky] HILC [full-sky] ILC [full-sky] 



1 


51.7% 


0(10"^) 


% 


3.7% 


2 


0(10-3)% 


0(10-3) 


% 


O(10-3)% 


3 


0.1% 


0.1% 




0.4% 


4 


0(10-3)% 


0(10-3) 


% 


0.1% 


5 


O(10-3)% 


0(10-3) 


% 


0.1% 


6 


O(10-3)% 


0(10-3) 


% 


%O(10-3) 


7 


1.2% 


0.1% 




1.3% 


8 


0.1% 


0.1% 




0.1% 


9 


O(10-3)% 


0(10-3) 


% 


0(10-3) 


10 


0.1% 


0(10-3) 


% 


0.1% 



TABLE IV: The probability (percentage) of occurrence of the 
multipole values 5*^ calculated from the data in the set {S'/} 
of values computed from MC Gaussian CMB maps. The data 
from the five-year full-sky foreground-reduced NILC, HILC 
and ILC maps were used. 

data were released, including a new version of the full-sky 
foreground- reduced ILC map jl^. We have considered 
this latest foreground-reduced ILC map, and performed a 
complete additional analysis of the Gaussianity of the five 
and seven-year versions of the ILC maps, whose details 
are given in Apenddix [Bj^ The main result of this ap- 
pendix is that the full-sky seven-year foreground-reduced 
ILC map also present a significant deviation from Gaus- 
sianity, which again is reduced substantially when the 
KQ75 mask is employed. In this way, our results are 
robust with respect to seven-year WMAP CMB data. 
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Appendix A 

Clearly from 1 000 MC maps one can calculate for 
each £ one thousand values of both S} and {i — 
1, • • • , 1 000 ) and the corresponding mean values Se and 
Ki. For the sake of brevity in what follows we focus on 
the skewness indicator S, but completely similar calcu- 
lations were used to have the probabilities for kurtosis 
indicator K. 



Note that there is no available seven-year HILC and NILC maps. 



£ NILC [KQ75] HILC [KQ75] ILC [KQ75] 



1 


76.7% 


34.1% 


73.7% 


2 


69.1% 


12.8% 


42.1% 


3 


4.6% 


3.7% 


5.7% 


4 


7.2% 


8.9% 


16.3% 


5 


2.7% 


2.2% 


11.3% 


6 


2.4% 


1.5% 


4.0% 


7 


51.4% 


57.6% 


54.1% 


8 


4.9% 


10.2% 


40.6% 


9 


7.0% 


13.8% 


55.0% 


10 


3.9% 


8.4% 


27.7% 



TABLE V: The probability (percentage) of occurrence of mul- 
tipole values Se calculated from the data in the set {S/} of 
values computed from MC Gaussian CMB maps. The data 
from the five-year NILC, HILC and ILC KQ75 masked maps 
were used. 



I NILC [full-sky] HILC [full-sky] ILC [fuU-sky] 



1 


0.6% 


0(10-3) 


% 


O(10-3)% 


2 


O(10-3)% 


0(10-3) 


% 


O(10-3)% 


3 


0.1% 


0(10-3) 


% 


0.1% 


4 


O(10-3)% 


0(10-3) 


% 


O(10-3)% 


5 


0.1% 


0(10-3) 


% 


O(10-3)% 


6 


0.1% 


0(10-3) 


% 


0.1% 


7 


0.4% 


0(10-3) 


% 


O(10-3)% 


8 


0.1% 


0(10-3) 


% 


O(10-3)% 


9 


0.1% 


0(10-3) 


% 


O(10-3)% 


10 


0.1% 


0(10-3) 


% 


O(10-3)% 



TABLE VI: The probability (percentage) of occurrence of the 
multipole values Ki calculated from the data in the set {-ft^/} 
of values computed from MC Gaussian CMB maps. The data 
from the five-year full-sky foreground-reduced NILC, HILC 
and ILC maps were utilized. 



i NILC [KQ75] HILC [KQ75] ILC [KQ75] 



1 


87.0% 


45.1% 


81.1% 


2 


2.3% 


1.9% 


23.6% 


3 


85.6% 


29.3% 


89.1% 


4 


2.3% 


9.0% 


57.1% 


5 


1.3% 


1.8% 


4.7% 


6 


82.6% 


4.2% 


98.1% 


7 


8.8% 


1.7% 


46.9% 


8 


96.6% 


93.2% 


86.6% 


9 


0.9% 


0.8% 


5.0% 


10 


65.0% 


54.2% 


97.7% 



TABLE VII: The probability (percentage) of occurrence of the 
multipole value Kt calculated from the data in the set {-ft'/} 
of values computed from MC Gaussian CMB maps. The data 
from the five-year foreground-reduced NILC, HILC and ILC 
KQ75 masked maps were employed. 



With the MC values and the mean Se one can cal- 
culate the percentages of values of the deviations \ S^ — Si\ 
calculated from 1 000 MC Gaussian maps which are 
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smaller than — Si\ with Si obtained from the data 
(full-sky and masked maps). For each multipole this 
number indicates how unlikely are the occurrences of the 
values obtained from the data (input maps) for that mul- 
tipole in the set of values {S*/} obtained from MC Gaus- 
sian maps. In this way one can calculate the probability 
of occurrence of a given multipole value Se (obtained from 
the data) in the set of MC values (obtained from the MC 
maps) for each foreground-reduced CMB map. In Ta- 
bles IIV| IVl IVII and IVIII we collect together the results of 
such calculations. 

Thus, for example, from Table HVl we have for the full- 
sky NILC, HILC, and ILC maps, respectively, the prob- 
ability of occurrence of the 5*6 values (in the set of MC 
values) is 0(10"^)%, whereas from Table IVll the proba- 
bility for Ka is, respectively, 0.1%, O{10-^)% and 0.1% 
for the full-sky NILC, HILC, and ILC input maps. 

The comparison of Table IW] with Table El and of Ta- 
ble |Vl] with Table Eni shows that the role of the KQ75 
mask is to cut down significantly the level of deviation 
from Gaussianity for all multipoles Si and Ki obtained 
from the foreground-reduced input maps. This is clear 
because the probabilities of occurrences for these multi- 
poles values in the set of MC multipole values increase 
substantially when the mask is employed. 

Although the estimates of probabilities collected in 
these tables give a clear quantitative indication of devia- 
tion from Gaussianity an overall assessment of the power 
spectra Se and Ki can be obtained through test of 
the goodness of fit for Si and Ki from the data as com- 
pared to the expected multipoles values obtained from 
the Gaussian MC maps. This point is discussed in Sec- 
tion ImBl 

Appendix B 

While we were in the final phase of writing a mod- 
ified version of this paper a new version of the full-sky 
foreground-reduced ILC map was released by the WMAP 
team . Since there is no available version of the NILC 
and HILC maps obtained from the seven-year WMAP 



data to be considered, here we present the results of a 
comparative analysis of deviation from Gaussianity per- 
formed by using the five and seven-year versions of the 
ILC as input maps. As the calculations are similar to 
those of Section IlII BI we refer the readers to that section 
for more details. 

Figure[5]shows the differential power spectra calculated 
from the full-sky five and seven-year foreground-reduced 
ILC input maps (ILC5 and ILC7, for short). Apart from 
some local deviation of the deviations \Si — Se\ and \Ke — 
Ki I this figure shows a deviation from Gaussianity, which 
is quantified in Table IVHll 

It is interesting to note that the deviation from Gaus- 
sianity as measured by our indicators is greater for the 
ILC7 than for the ILC5 input map. Concerning this point 
some words of clarification are in order here. First, we 
note that the details of the algorithm used to compute 
the ILC7 maps are the same as those of the ILC5 map. 
However, to take into account the most recent updates to 
the calibration and beams, the frequency weights for each 
of the 12 regions (in which the sky is subdivided in the 
ILC method) are slightly different in the calculation of 
the ILC7 map. Second, the difference between the ILC7 
and ILC5 maps is a map whose small-scale differences 
are consistent with the pixel noise, but with a large-scale 
dipolar component, with the large-scale differences being 
consistent with a change in dipolc of 6.7 /iK[l5i. Thus, 
the resultant ILC 7 map is not indistinguishable from the 
ILC5 map, and the differences between them have been 
captured by our indicators. 

Figure[S]shows the differential power spectra calculated 
from a five-year and seven year version of the foreground- 
reduced ILC maps with a KQ75 mask. This figure along 
with Fig. [5] show a significant reduction in the level of de- 
viation from Gaussianity when both ILC5 and ILC7 are 
masked. To quantify this reduction we have calculated 
X^/dof for these input maps with the KQ75 mask, and 
have collected the results in Table IIXI The comparison 
of Table IVIIII and Table IIXI shows quantitatively the re- 
duction of the level of Gaussianity for the case of CMB 
masked maps. 
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FIG. 5: Differential power spectrum of skewness \Si — Si\ (left) and kurtosis \Ki — Ki\ (right) indicators calculated from 
tlie full-sky foreground-reduced ILC input map obtained from the WMAP five-year and seven-year data. The 68% and 95% 
confidence levels are indicated, respectively, by the dashed and dash-dotted lines. 





for St [full-sky] 


for Ki [full-sky] 


ILC5 


35.7 


2368 


ILC7 


103 


10507 



TABLE VIII: Results of the test to determine the goodness of fit for St and Ki multipole values calculated from the 
foreground-reduced ILC5 and ILC7 full-sky maps as compared to the expected multipoles values from the MC Gaussian full- 
sky maps. 



Itf- 



1 0" 



1 0"" 



1 0" 



S„ : ILC — Syr, ILC — 7yr KQ7S masked maps 



1 0" 



O = ILC-7yr 
• = ILC-5yr 



1 2 3 



4 5 6 7 
multipoles (t) 



9 1 



1 c- 



1 0" 



1 0" 



1 0" 



10" 



K„ : ILC — Syr, ILC — 7yr KQ7S masked mops 



^ = ILC-7yr 
• = ILC-5yr 



^ ^ -. 



2 3 4 5 6 7 
multipoles (0 



8 9 10 



FIG. 6: Differential power spectrum of skewness \Si — Se\ (left) and kurtosis \Ki — Ke\ (right) indicators calculated from the 
five-year and seven-year foreground-reduced ILC input maps with a KQ75 mask. The 68% and 95% confidence levels are 
indicated, respectively, by the dashed and dash-dotted lines. 





X" for Se [KQ75] 


X" for Ke [KQ75] 


ILC5 


1.2 


0.4 


ILC7 


0.7 


0.2 



TABLE IX: Results of the x^ test to determine the goodness of fit for Si and Ki multipole values calculated from the foreground- 
reduced ILC5 and ILC7 maps with a KQ75 mask as compared to the expected multipoles values from the Gaussian MC masked 
maps. 
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